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We investigate the femtosecond kinetics of an optically excited 2D magneto-plasma at interme- 
diate and high densities under a strong magnetic field perpendicular to the quantum well (QW). 
We assume an additional weak lateral confinement which lifts the degeneracy of the Landau levels 
partially. We calculate the femtosecond dephasing and relaxation kinetics of the laser pulse excited 
magneto-plasma due to bare Coulomb potential scattering, because screening is under these con- 
ditions of minor importance. In particular the time-resolved and time-integrated four-wave mixing 
(FWM) signals are calculated by taking into account three Landau subbands in both the valance 
and the conduction band assuming an electron-hole symmetry. The FWM signals exhibit quan- 
tum beats mainly with twice the cyclotron frequency. Contrary to general expectations, we find 
no pronounced slowing down of the dephasing with increasing magnetic field. On the contrary, one 
obtains a decreasing dephasing time because of the increase of the Coulomb matrix elements and the 
number of states in a given Landau subband. In the situation when the loss of scattering channels 
exceeds these increasing effects, one gets a slight increase at the dephasing time. However, details of 
the strongly modulated scattering kinetics depend sensitively on the detuning, the plasma density, 
and the spectral pulse width relative to the cyclotron frequency. 



42.50.Md, 42.65.Re, 42.50.-P, 78.47.+p 
I. INTRODUCTION 

Femtosecond pulse excitation in semiconductors in- 
duces transient carrier populations, which can be stud- 
ied through ultrashort-pulse nonlinear optics to elucidate 
many-body effects, such as time-dependent Coulomb cor- 
relations. Numerous experimental and theoretical stud- 
ies have been devoted to this problem-in the magnetic- 
field-free case in the last ten yearslHj Studies of fem- 
tosecond laser spectroscopy in the presence of a strong 
magnetic field are relatively rare and are mainly focused 
on magneto-excitons at low excitation densitiesErHO. 
For a strong resonant laser pulse which excites an 
electron-hole (e-h) system with a density above the Mott 
ionization density, the kinetics is dominated by Coulomb 
scattering in the correlated e-h plasma or in the cor- 
related magneto-plasma if a strong magnetic field is 
present. The femtosecond quantum kinetics and the ex- 
pected FWM signal of the non-equilibrium e-h plasma 
have been treated recently by Vu et aln using bare 
Coulomb potential for times shorter than the build- 
up time of screening. While the density-dependence of 
the optical spectra of a quasi-equilibrium 2D magneto- 
plasma have been calculated,El the kinetics of a magneto- 
plasma has been studied neither experimentally nor the- 
oretically. Moreover, the existing theory for the kinetics 
with a magnetic field is also not yet developed as far as 
in the field-free case. However, as experimental studies 
of the relaxation and dephasing kinetics in QW's and 
superlattices are in progresalJ, we will provide a first 
calculation of the femtosecond FWM signal of a reso- 
nantly excited magneto-plasma in a single quantum well. 
We concentrate on the study of carrier-carrier scattering, 



while phonon and disorder scattering is not considered. 

A strong magnetic field perpendicular to the QW plane 
forces the carriers on cyclotron orbits, thus causing an 
additional quantum confinement. In a magnetic field of 
the order of more than 10T only a few Landau levels 
both in the conduction and valence band have to be con- 
sidered. The high degeneracy of these Landau levels is 
partially lifted by spatial inhomogeneities and by size ef- 
fects. The numerical treatment of the broadening due 
to disorder, e.g. interface fluctuations, would require a 
stochastic averaging over many FWM signal calculations 
which is beyond present day numerical possibilities for 
the complex carrier kinetics treated here. Therefore, we 
lift the degeneracy partially by a weak parabolic confine- 
ment potential. The Landau levels are broadened by this 
confinement into Landau subbands with a band width 
of about 2 meV. This weak parabolic confinement allows 
to treat the single-particle Schrodinger equation still ex- 
actly in the presence of the magnetic fieldO These eigen- 
functions of the 2D electron provide a convenient expan- 
sion basis for the non-equilibrium many-body problem. 
A strong quantum confinement — here caused by the 
QW and the strong magnetic field — is generally be- 
lieved to slow down the relaxation kinetics because of tha 
reduction of the phase space for scattering processes.E£l 
Indeed, for low-density magneto-excitons an increase in 
dephasing time has been observed.E3 However, this rule 
does not apply to Coulomb scattering, as has been shown 
recently for the example of exciton-exciton scattering 
in QW wires. til For decreasing wire width a reduction 
of the dephasing time has been found and explained in 
terms of the increase of the Coulombic interaction ma- 
trix elements which overcompensated the reduction of 
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phase space for the scattering processes. Therefore it is 
not obvious, how the Coulomb intra- and inter-Landau- 
subband scattering will influence the resulting dephasing 
time in a dense magneto-plasma in detail. 

We present a kinetic study for a femtosecond laser- 
pulse excited 2D dense non-equilibrium magneto-plasma 
in a QW in the framework of the semiconductor Bloch 
equations combined with Coulomb scattering rates (Sec. 
II). We expand the density matrix of a two-band semi- 
conductor in the eigenfunctions of the 2D electron in 
the presence of the strong magnetic field and the weak 
parabolic confinement. We formulate the scattering 
terms for the population distribution functions of the 
various Landau subbands and for the optically induced 
polarization components between the Landau-subbands 
in the conduction and valence band in the form ol non- 
Markovian quantum kinetic scattering integralso and 
in the form of semiclassical Boltzmann-type scattering 
rates. Pronounced quantum kinetic effects are expected 
for time scales shorter than typical inverse frequencies. 
For the considered high magnetic fields, the cyclotron fre- 
quencies are so large that quantum kinetic effects should 
be of minor importance. Furthermore, we prefer the use 
of the simpler semiclassical kinetics, because there exist 
still some conceptual difficulties for the quantum kinetic 
description, connected with the spectral electron Green 
functions in a strong magnetic field, as will be discussed 
below. We calculate the time-resolved (TR) and time- 
integrated (TI) four-wave mixing (FWM) signals for two 
50 fs pulses by taking into account up to three Landau 
subbands in both the valance band and the conduction 
band. The carrier frequency of the two delayed pulses 
is tuned slightly above the un-renormalized energy gap. 
We simplify the problem by assuming equal effective elec- 
tron and hole masses, as it can be approximately realized 
in strained QWs. Naturally, unequally effective masses 
will lead to more complicated quantum beat structures in 
the FWM signals and modify to some extend also the re- 
sulting relaxation and dephasing rates. Thus our present 
study should be seen only as a first idealized model calcu- 
lation. The detailed numerical results based on the semi- 
conductor Bloch equations with Boltzmann-type scatter- 
ing kinetics are presented in Sec. Ill for intermediate and 
high plasma densities. The FWM signal is calculated by 
an adiabatic projection technique which is appropriate 
for thin samples where propagation effects are not im- 
portant. For many conditions pronounced quantum beat 
structures with two times the cyclotron frequency are 
obtained both in the TR and TI FWM signals. Further- 
more, often relaxation oscillations are seen between the 
populations of two subbands, before they relax to their 
stationary values. The dephasing times obtained from 
the TI FWM signal are surprisingly short and strongly 
modulated as a function of the magnetic field. The un- 
expectedly short resulting dephasing times between 100 
and 200 fs depend not only on the field values, but also 
on the detuning, the pulse width and the excited carrier 
densities. In spite of the remaining conceptual difficul- 



ties we present also first quantum kinetic calculations 
and compare the FWM signals with those obtained with 
Boltzmann-type kinetics in Sec. IV. A conclusion of our 
main results is given in Sec. V. In an appendix we give 
a list of the analytically calculated Coulomb matrix ele- 
ments for the first three subbands. 



II. MODEL AND KINETIC EQUATIONS 

A. Model and Hamiltonian 

We start our investigation of a quantum well in the 
x-y plane which is further restricted by an additional 
weak lateral confinement in x direction which lifts the 
degeneracy of the Landau levels partially and provides a 
weak inhomogeneous broadening. The confinement is as- 
sumed to be given by a harmonic oscillator potential for 
the lower-lying states. This model has been used with 
a strong additional confinement for the theory of the 
optical-properties of thermal magneto-plasmas in QW 
wiresJli A strong magnetic field B > 10T is applied per- 
pendicular to the well. We assume electron-hole sym- 
metry, i.e. m e = rrih = m with m e (rrih) denoting the 
effective electron (hole)mass. Within the Landau gauge 
A = xBe y and effective mass approximationEj to the 
lowest QW subband, the stationary 2D Schrodinger equa- 
tion of a single electron can be solved exactly (U = c = 1): 
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with H n (x) standing for the n-th order Hermite polyno- 
mial and a = \/mQ x is the inverse of the amplitude of 
the zero-point fluctuations. Q x is the effective oscillation 
frequency fl x = y^fl 2 + uj 2 , where oj c = eB/m is the cy- 
clotron frequency and f2 is the frequency of the additional 
confinement potential. The shift §Xf- — —u> c k/(m£l 2 ) re- 
sults from the balance between the Lorentz force and the 
harmonic restoring force e(k/m)B = mtt 2 Sxk. 

The single-particle energy spectrum of a conduction 
band electron in n-th Landau level and with momentum 
k is given by E cnk = E g /2 + e nk with 
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and Eg standing for the energy gap. The dependence of 
the energy on the momentum in y-direction results from 
the additional confinement of the QW, which broadens 
each Landau level into a small Landau subband. The 
momenta are restricted to 



|fc| < L x mfi*/(2w c ) 



(5) 



because the centec of cyclotron should lie within the 
sample width L x Jl3 It is interesting to notice that when 
lu c 3> the inhomogeneous broadening is of the same 
for all magnetic fields as the increase of k space in Eq. 
(0) is totally compensated by the increase of the effective 
mass in Eq. (^). 

With Coulomb interaction and with the interaction 
with a coherent classical light field, the many-body 
Hamiltonian for the electrons in the conduction and va- 
lence bands is in the basis of above given magnetic eigen- 
functions 

H = ^ E unk c\ nk C in k + i ^2 Vni;jm (Q, k, k') 

vnk v\v-2\nmiy.kk' q 

X ^vtnk+q^ik'-qC^jk'C-vxmk + Hj , (6) 

with v = c and v standing for the conduction band and 
the valence band respectively. E vnk — —E cn k is the en- 
ergy spectrum of an electron in the valence band. 

The Coulomb interaction matrix elements 
V n i;jm(Q>k)k f ) describes the scattering of an electron 
with momentum k from m-th Landau subband to n-th 
subband with momentum k + q and an electron with mo- 
mentum k! from j-th subband to «-th with momentum 
k' — q. Due to the assumed e-h symmetry, these four 
Landau subbands can be either in the conduction band 
or in the valence band. The matrix elements are given 
in terms of the shifted magnetic eigenfunctions by the 
following integral 
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e is the background dielectric constant. After substitut- 
ing x — ► x + Sxk and x' — ► x' + 6xk> in Eq. (Q), one 
gets 
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+ A(fc - k')]}<p* n (x + Sx q )ct)*(x' - 6x q )(j> j (af)4> m (x) , (8) 

with A = —muj c /a 2 . After an approximation concern- 
ing the momentum dependence, the integrals of Eq. (0) 



may be carried out analytically. In the appendix the cor- 
responding approximation and the final expressions of 
V n i:j m (q,k,k r ) in terms of zeroth- and first-order modi- 
fied Bessel functions are given. 

Hj in Eq. (^) denotes the dipole coupling with the light 
field E(t). In the assumed e-h symmetry it contains only 
transitions between Landau subbands of the same order: 
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In this equation, d denotes the optical-dipole matrix ele- 
ment. The light field is further split into E(t) = E (t)e luJt 
with to being the central frequency of the coherent pulse. 
Eo(t) describes a Gaussian pulse E^e - * ^ St with 5t de- 
noting the pulse width. 



B. Kinetic Equations 



Following the same scheme as for the magnetic-ficld- 
free case in Ref. ||, appropriately modified, one may 
build the semiconductor Bloch equations for the QW in 
a strong perpendicular magnetic field B in the basis of 
the magnetic eigenfunctions as follows: 

Pv^n.v' ,n' ,k — Pu,n,u' ,n' ,k Icoh P v ^ n ^ v ' : n ' I SCatt ' (^^) 

Here p v . n y ,n' ,k represents the single-particle density ma- 
trix. The diagonal elements describe the carrier distri- 
bution functions p u ,n,v,n,k — funk of the n-th Landau 
subband and the wavevector k as diagonal elements, and 
the off-diagonal elements describe the interband polar- 
ization components, e.g. p c ,n,v.n,k — Pnke~ lult - For the 
assumed e-h symmetry, f enk = f hnk = f nk and the po- 
larization has only components between subbands of the 
same quantum number n in the conduction and valence 
band, which simplifies the problem considerably. 

The coherent part of the equation of motion for the dis- 
tribution function is in the rotating wave approximation 
given by 
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The first term describes the generation rate by the laser 
pulse, d is the optical dipole matrix clement. For the as- 
sumed e-h symmetry transitions between different sub- 
band quantum numbers are not allowed. The second 
term describes the exchange interaction correction of the 
exciting laser by the e-h attraction, thus it can be seen as 
a local field correction of the time-dependent bare Rabi 
frequency dEo(t). The retarded quantum kinetic scat- 
tering rates for the coH&idered bare Coulomb potential 
scattering are given bytrB 



3 



scatt 
t 



-8 \ V m;j m (<l,k,k' + q)f 



qk'mij 



X 



f dt'Re{e l - l{e ™ k -i- E '> k+E i*' 

J — OO 



+q -s ik ')-r](t-t') 



X [fnk(t')fik'(t')(l - f m k-S'W - fjk'+ q (t')) 

- f,nk- q (t')f jk ' +q (t')(l ~ fnk(t'))(l ~ fik'(t')) 

- P; k {t')P mk . q {t'){f ik ,{t') - fjk- +q {t')) 

- P jk , +q {t')P* k ,{t'){f nk {t') - f mk - q (t'))]} . (12) 

In the derivation of this formula vertex corrections have 
been neglected. The scattering self-energy is evaluated 
in time-dependent RPA. Successively, the Coulomb po- 
tential is taken as a bare instantaneous potential and 
the two-time-dependent particle propagators have been 
expressed in terms of single-time density matrix ele- 
ments and retarded (or advanced) Green functions us- 
ing the generalized Kadanoff-Baym ansatz. For simplic- 
ity the retarded (and advanced) Green functions have 
been approximated by a diagonal free-particle Wigner- 
Weiskopf form, i.e. G r nk (t,t') = G a n * k {t' ,t) = -i®{t - 



-i(E g /2+e nk )-~f](t-t') 



The effective damping T in Eq. 



( |12| ) is the sum of the four imaginary parts of the re- 
tarded electron self-energy which is assumed to be sim- 
ply the damping constant 7, so that T ~ 47. We will 
discuss later that this approximation for the imaginary 
part of the retarded e-self-energy leads to some unphys- 
ical features in the quantum kinetics of the considered 
magneto-plasma. For further progress in the quantum ki- 
netics of a magneto-plasma, the calculations of the spec- 
tral functions have to be made more self-consistent in 
order to include the retarded onset and the magnetic 
field dependence of damping and the band mixing by the 
coherent light pulses and the mean-field Coulomb inter- 
actions. These refinements have been developed already 
for the simpler case of n scattjering with optical phonons 
without magnetic fieldHLStj However, for the time be- 
ing the complexity of the Coulomb quantum kinetics in 
a magneto-plasma prevents us to include these improve- 
ments in the present numerical evaluation. 

The coherent time evolution of the interband polariza- 
tion components are 
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The first term gives the free evolution of the polarization 
components with the detuning 

S n (k) =2s n (k)- A - 2 J2Vnm;nm(q,k,k)f mq {t) (14) 

with Ao = u) — Eg . Ao is the detuning of the center 
frequency of the light pulses with respect to the unrenor- 
malized band gap. The second term in Eq. (O describes 



again the excitonic correlations in the magneto-plasma, 
while the final factor describes the Pauli blocking. 

The dephasing of the polarization components are de- 
termined by the following quantum kinetic scattering in- 
tegral 
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N n ji(k, k', q, t') = (1 - /„ fc (0)(l - fiki +q {t'))fik- 
+ fnk(t')(l - fik'{t'))fjk' +q {t') - P jk > +q {t')P* kl {t') , (16) 

and the energy difference S n (k) = 2e n (k) — Aq — 

^ Em,o Km;nm(9i k, k)f mq (t). Aq = 



Eg is the detun- 



ing with respect to the unrenormalized band gap. Here 
T2 is introduced phenomenologically to describe addi- 
tional slower scattering processes. Eqs. @-(|l|), (|l|) 
and (|l5|) are the quantum kinetic Bloch equations for a 
magneto-plasma. 

In the longtime limit the quantum kinetic scattering in- 
tegrals be transformed into Markovian Boltzmann-type 
scattering rates by pulling the slowly varying distribu- 
tions and polarization components outside the scattering 
integrals at the upper limit t and by replacing the mem- 
ory kernels / ^ dt' exp{ [-i{e mk -q - e nk +e jk , + q - e ik < ) - 
T](t-t')} in Eq. @ by 2n5{e mk - q - e nk + £jk'+ q - £%k>) 
and J ^ dt' exp{[-i(£ mfe _,j + e nk + e jk >+ q - e iV - A ) - 
T](t~t')} in Eq. |[|) by 2TTS(e mk - q + e nk + e jk ' +q -e lk ' - 
A ) in the limit of vanishing damping. In our calcula- 
tions we use for numerical convenience energy resonances 
with a small Gaussian damping, which does not lead to 
such severe violations of the energy conservation as the 
Lorentzian resonances do. 

It is noted that besides the e-h symmetry, we have fur- 
ther assumed that Coulomb Auger scattering across the 
gap can be neglected. 



III. NUMERICAL RESULTS BASED ON 
BOLTZMANN KINETICS 

We perform a numerical study of the Bloch equations 
in the Boltzmann limit to study TR and TI FWM signals 
for high magnetic fields (B > 10 T) for two degenerate 
Gaussian pulses of a width of 50 fs and a variable de- 
lay time t. The pulses travel in the direction ki and 
k2. We use an adiabatic projection technique in order 
to calculate the polarization in the FWM direction with 
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wavevector 2k2 — ki described in detail in Ref. 16, This 
technique is suitable for optically thin crystals, where the 
spatial dependence can be treated adiabaticallytZl To do 
so, we replace the single-pulse envelope function in Eq. 
(y) by two delayed pulses E (t) = E {t) + E (t - r)e lv 
with the relative phase ip — (k.2 — ki) • x resulting from 
the different propagation directions. The projection tech- 
nique is used with respect to this phase. We further as- 
sume a resonant excitation with a fixed excess energy of 
Ao = 26.4 meV. The intensity of each pulse will be given 
by dE (t)dt — x 7 ?, where x denotes the fraction of 
a 7r-pulse defined without local field corrections. Except 
to /i which is assumed to be equal to m e , all the other 
material parameters are taken for bulk GaAs with an ex- 
citon Rydberg of 4.2 meV and a Bohr radius of 14 nm. 
For the sake of weak confinement, the harmonic oscilla- 
tor frequency fl is assumed to be 0.1u> c for B = 7 T, 
corresponding to Q = 1.2 meV. By fitting the energy f2 
with the energy splitting to the two lowest subbands of 
a rectangular square well 3tt 2 / (2mL 2 ) , we get the width 
in the additional confinement direction to be L T = 118 
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FIG. 1. The pulse intensity spectrum \E(ui)\' 2 plotted in 
both linear and log scale for a 50 fs pulse, together with the 
unrenormalized energybands for the optical transitions be- 
tween Landau subbands n(—0, 1, 2) plotted as solid lines for 
B = 12 T and as dashed lines for B = 18 T. 



In Fig. 1 we plot the pulse spectrum together with the 
three transitions bands for two magnetic fields of B = 12 
T (solid lines) and 18 T (dashed lines). The first transi- 
tion band corresponding to Pq is due to transitions be- 
tween the lowest Landau subband n = in the valance 
band to the lowest subband n = in the conduction 
band. This transition band starts an effective frequency 
£l x above the band gap at least in the low-density limit. 
Similarly the second transition band for the n=l Landau 
subbands in both valence and conduction bands requires 






FIG. 2. The distribution functions of three Landau sub- 
bands /o, /i and /2 are plotted as functions of t and k for a 
one-pulse excitation with B — 12 T. The units of k in the fig- 
ures are one tenth of the maximum allowed value of k defined 
in Eq. (§). 



the excess energy 3f2 x . The Hartree-Fock terms in Eq. 
( |l4| ) may red-shift the energy of transition for the ini- 
tial stage of excitation. From Fig. 1 we see that for high 
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magnetic fields, the pulse mainly populates in the lowest 
Landau subband. 
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FIG. 3. Electron densities of three Landau subbands Ni 
and the total density Nt and the incoherently summed polar- 
ization plotted against time t for B = 12 T (Fig. 3(a)) and 18 
T (Fig. 3(b)). 
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FIG. 4. TR and TI FWM signals versus time t and delay 
r, respectively, for B = 12 T (Fig. 4(a)) and 18 T (Fig. 4(b)). 
For TR FWM, the solid curve is for r = 120 fs and dashed 
curve for 240 fs. The dotted line is the exponential fit to the 
TI FWM signal. Note that the scale of the delay time r for 
TI FWM signals (top frame) is different from the time scale 
of the TR FWM signals (bottom frame). 



As mentioned before, we include in our calculations 
three Landau subbands in the valence and conduction 
band. It is important to include higher Landau subbands 
for the kinetics after the excitation. In principle it is true 
even in situations where the pulses mainly excite e's and 
h's only in the pair of the lowest Landau n=0 subbands: 
A certain number of e's excited in the n=l Landau sub- 
band can only relax to the n=0 subband, if simultane- 
ously the same number of e's is excited from the n=l 
to the n=2 subband if only three subbands are included. 
On the other hand, say if one includes up to 5 Landau 
subbands, then it is possible to only scatter one electron 
from the second Landau subband to the fifth subband 



and in the mean time to release three electrons from 
the second subband to the lowest subband. Therefore, 
a too small number of Landau subbands may change the 
Coulomb kinetics. For this reason, we restrict ourselves 
to the high magnetic field regime in which for the cho- 
sen detuning e's and h's are excited mainly in the lowest 
n=0 Landau subbands. More Landau subbands are nec- 
essary in order to extend the kinetics to lower magnetic 
fields or for larger detunings. However, the expansion of 
number of matrix elements V n i-j m increases as iV 4 with 
N being total number of Landau subbands considered. 
With 7Y = 3 in our model, the number of form factors is 
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already 81, while it will be 256 with N = 4. 



A. Intermediate-density case 

We first discuss the Boltzmann kinetics for an inter- 
mediate excitation density. To do so, we choose a pulse 
with x = 0.1. We first show the distribution functions 
fnk(t) for n = 0, 1 and 2 versus t and k for a one-pulse 
excitation in Fig. 2 for B = 12 T. Here and hereafter, 
the additional transverse relaxation time T-x is taken as 
300 fs. The energy conserving (^-functions in the col- 
lision terms are replaced by Gaussian functions with a 
width a = 0.66 meV which is much smaller than both io c 
and the inhomogeneous broadening. Note again that k 
is limited by Eq. (|5|). One sees how the carriers relax at 
later times in the lower fc-states of the Landau subbands. 
We plot the various carrier densities Ni(t) = ^2 k fik(t) 
against time t for the magnetic fields B = 12 T (Fig. 
3(a)) and 18 T (Fig. 3(b)). For B = 12 T, one sees 
pronounced relaxation oscillations of the populations of 
the lowest subbands No(t) and Ni(t) in the first 400 fs. 
The period of these relaxation oscillations (w 60 fs) with 
successive overshoots in the population of the Landau 
subband n = and n = 1 is given by the characteris- 
tic interband Coulomb scattering rate between these two 
lowest subbands. The total carrier density N t = J2i Ni 
shows a slight Rabi flopping overshoot. The incoherently 
summed polarizations P(t) = ^2 ik \Pik \ for the two fields 
are also plotted in the same figures. Strong quantum 
beats with frequency 2lo c are revealed in the polarization 
for lower magnetic fields up to 16 T. This effect may be 
measured by FWM. For B > 16 T, both the population 
oscillation and the beating of the polarization amplitude 
die out (see Fig. 3(b)). Moreover, one sees that in the 
latter case the carriers are excited practically only in the 
lowest Landau subband. 

In Fig. 4 we plot typical TR FWM signals versus time 
t for B = 12 T (Fig. 4(a)) and 18 T (Fig. 4(b)) for two 
pulses with the delay times r = 120 fs (solid curve) and 
240 fs (dashed curve) , respectively. We find strong beat- 
ing for lower magnetic fields, whereas for higher fields 
the beating becomes weaker. Only very small beats are 
found for B = 18 T in Fig. 4(b). The beating totally 
disappears for B = 20 T. The beating frequency for the 
TR FWM signal is exactly 2lo c for all magnetic fields 
where beating exists. From the TI FWM signal as a 
function of the delay time r one can obtain the effective 
dephasing time. Therefore we plot the TI FWM signals 
in the same figures, however, projected on the right and 
top scales. Again one observes quantum beats in the TI 
FWM signals for lower fields. Similar beating effect has 
been obtained based on a three-level model in explaining 
heavy hole- light hole beats by Leo et al.&3 Moreover, the 
beats become weaker when the magnetic field increases 
because then the second subband becomes less impor- 
tant. For B = 16 T the beating has already disappeared 



in TI FWM which means that the polarization compo- 
nents P\k are negligible in comparison to Pok ■ This 
is consistent with Fig. 3(b), where the electron (hole) 
population in the second Landau subband is extremely 
low and the beating structure disappears in the incoher- 
ently summed polarization. The beating frequency for 
the lower field is also 2ui c . Note again the different scales 
used for t and r. The dotted lines in both figures are 
the exponential fits to the TI FWM signal. From the 
fit one gets the effective dephasing time which is plotted 
as a function of the magnetic field B in Fig. 5. Interest- 
ingly, we find for large magnetic fields that the dephasing 
time decreases with increasing field forlOT<_B<16T. 
However, there is a transition at 16 T. For B > 16 T the 
dephasing time increases with increasing field, although 
much slower than it decreased before. 
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FIG. 5. Dephasing time as a function of B. 

For a fixed pulse, several effects compete with each 
other when the magnetic field increases. On one hand, 
the number of Landau subbands which contribute to the 
Coulomb scattering kinetics decreases. In particular the 
contributions to the dephasing from the intra- and inter- 
subband scattering of the higher Landau subbands as 
well as the inter-subband scattering between the higher 
and lower subbands decrease. For large populations in 
one subband the Pauli blocking may further reduce also 
the intra-band scattering rates. All these effects increase 
the dephasing time. On the other hand, with increas- 
ing B field the degeneracy of Landau subbands increases 
and the matrix elements of the Coulomb scattering of 
Eq. (^) become larger. How an increasing degeneracy 
increases the scattering rates can be seen from Eqs. ( |l2] ) 
and ( |l5| ) in the extreme limit where the confinement is 
lifted. Then there is no k and k' dependence in either 
distribution functions or polarization functions. There- 
fore X)fc'rjf replaced by Smuj c / (2ir) (S standing for the 
2D area)li3 which increases as B increases. Such an ef- 
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FIG. 6. Distribution functions of three Landau subbands 
fo, fi and f2 versus t and k for one-pulse excitation for B = 12 
T. The units of k in these figures are one tenth of the maxi- 
mum allowed value of k defined in Eq. ^ . 



increase in the scattering rates is overcompensated by the 
loss of the contributions of the inter-subband scattering 
processes and the reduction of the intra subband scatter- 
ing due to Pauli blocking, so that a net slow increase of 
the dephasing time with increasing B field is obtained. 
For 10T < B < 16T the increasing degeneracy and the 
increasing Coulomb matrix elements dominate over the 
reduction of inter-subband scattering so that the dephas- 
ing time decreases with increasing B field. 
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FIG. 7. Electron density and incoherently summed polar- 
ization are plotted against time t for B = 12 T (Fig. 3(a)) 
and 18 T (Fig. 3(b)). N n (n = 0, 1, and 2) is carrier density 
in n-th Landau subband. N t = Ni + N2 + N3 is the total 
density. 



feet is also kept after the weak broadening in k space is 
included. Thus both the increased degeneracy and the 
increased Coulomb matrix elements reduce the dephas- 
ing time. Particularly, for B > 16T the above discussed 



B. High-density case 

In order to obtain a dense magneto-plasma with strong 
carrier interaction, we choose \ — 0.3 corresponding 
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roughly to a 7r/3 pulse, which provides a very high ex- 
citation density. The detuning is kept to be Ao = 26.4 
meV. Note that the local field renormalization will ac- 
tually turn the optical Bloch vector through an angle 
considerably larger than ir/3. We show the distribution 
function f n k{t) with n = 0, 1 and 2 versus t and k for 
one-pulse excitation in Fig. 6 for B = 12 T. 



r(fs) 



100 120 140 160 180 200 220 240 260 280 



■ III 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 


■ | i i i i | i i i i | i i i i | i i i i 






B = 12T 




'i V - . 
A' \ 






1 11 






l > < 








(a) > 




i y \/^ 

(..!.,....,.. 





3 

10 4 <^ 



200 400 600 800 1000 1200 1400 
t(fs) 

r(fs) 

100 120 140 160 180 200 220 240 260 280 



h 




600 800 1000 1200 
t(fs) 

FIG. 8. TR and TI FWM signals versus time t and delay 
t respectively for B = 12 T (Fig. 4(a)) and 18 T (Fig. 4(b)). 
For TR FWM, the solid curve is for r = 120 fs and dashed 
curve for 240 fs. The dotted line along TI FWM signal is the 
exponential fit to the TI FWM signal. 

From the figure one can see that the carrier distribu- 
tion function for the lowest Landau subband is close to 
0.5. This number approaches 1 when B is around 18 T. 
One sees a pronounced Rabi-flopping in the two lowest 
Landau subbands. We further plot the carrier density 
against time t for two different magnetic fields B = 12 
T (Fig. 7(a)) and 18 T (Fig. 7(b)). The incoherently 



summed polarizations for the two fields are also plot- 
ted in the same figures. Differing from the intermediate- 
density case, strong quantum beats are revealed for both 
field strengths, while the population relaxation oscilla- 
tions with a period of about 100 fs arc again only well 
resolved for the low-magnetic field case (B = 12 T). 
Again, the lowest Landau subband is always excited most 
strongly. 

In Fig. 8 we plot again typical TR FWM signals against 
time t for B = 12 T (Fig. 4(a)) and 18 T (Fig. 4(b)) 
for two pulses with the delay times r = 120 fs (solid 
curve) and 240 fs (dashed curve), respectively. Similar 
to the intermediate-density case, we find quantum beat- 
ing for fields as high as 20 T with a frequency of ex- 
actly 2u> c , although for large fields the beating become 
weaker. The TI FWM signals are plotted in the same 
figures as functions of delay time t (upper time scale). 
Again one observes quantum beats in the TI FWM sig- 
nals. Interestingly, the beating frequency is close to that 
in the incoherently summed polarization and is a little 
bit smaller than 2co c for higher fields. For example, when 
B = 18 T, it/uj c = 66.5 fs whereas from the beating in 
TI FWM, one finds a period around 76 fs. This differ- 
ence may be understood in terms of many-body effects, 
such as band gap renormalization, Pauli blocking and 
excitonic enhancement. We conclude that the quantum 
beats in the high-density case are also caused by interfer- 
ences between the optically induced polarizations P n =ok 
and P n= i.k of the two lowest subbands. 

The effective dephasing is plotted as a function of the 
magnetic field B in Fig. 9. Once more we find a grad- 
ual overall decrease of the dephasing time with increasing 
field in the regime from 10 to 20 T. This is in agreement 
with our earlier results in the intermediate-density case 
for 10 T< B < 16 T. The dephasing is mainly controlled 
in the whole range by the kinetics in the lowest two Lan- 
dau subbands. Compared to the intermediate-density 
regime, one observes a sharp dip around 15 T in the de- 
phasing time. By investigating the distribution function, 
we find the dip around 15 T is mainly due to the function 
of the second pulse. When B is smaller than 15 T, the 
second pulse further pumps electrons from valence band 
to conduction band. Due to the density dependent red 
shift of the subbands and the Pauli blocking for the op- 
tical transition in the lowest subband, the second pulse 
increases the population of the second subband relative to 
that of the first one. The ratio of the densities in the low- 
est two subbands N /Ni falls from around 5.7 after the 
first pulse to 4 after the second pulse when B=14 T, and 
that number goes from 8 to 3.9 for B = 15 T. This means 
the role of second subband becomes more important to 
the dephasing. This gives rise to the sharp decrease of 
the dephasing time. However for B > 15 T, the lowest 
Landau subband is already highly populated and due to 
the high magnetic field value, the second Landau sub- 
band is already relatively far away from the center of the 
pulse, so that the second pulse depopulates the electrons 
from the first conduction subband. Although this makes 
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again the density of lowest subband closer to the second 
one, the ratio N /Ni after the second pulse goes from 3.9 
at 15 T to 6 at 16 T. This clearly indicates the loss of the 
second subband to the dcphasing processes and explains 
the fast increase of the dephasing time from 15 to 16 T. 
After 16 T, although the first pulse populates more carri- 
ers to the lowest subband with increasing magnetic field 
due to the density-dependent band gap shift, the rate of 
change becomes smaller as the lowest subband is nearly 
filled. However, the depopulation due to the second pulse 
becomes stronger with increasing magnetic fields. This 
can be seen from the fact that No/N\ after the second 
pulse goes down from 6 for B = 15 T to finally 5 for 
B = 20 T. Our results show that the decrease of Pauli 
blocking together with the increasing of intra-subband 
scattering matrix element in the lowest subband domi- 
nate in the high magnetic field range the dephasing and 
explain the resulting decrease of dephasing time. 




14 16 
B(T) 

FIG. 9. Dephasing time as a function of B. 



IV. COMPARISON WITH QUANTUM KINETIC 
CALCULATION 



limit and we find that the results, particularly for the 
TI FWM signals and therefore the dephasing are highly 
sensitive to the actual value of the damping constant. 

In Fig. 10 we plot the TR FWM signal for pulses with 
X = 0.3 and r = 240 fs for B = 18 T. The solid curve 
is based on quantum kinetics with T = 1.32 meV. The 
dashed curve is the prediction based on Boltzmann ki- 
netics as shown in Fig. 8(b). From Fig. 10 one can see 
that the two theories give qualitatively similar TR FWM 
signals, although the coherence of the polarization com- 
ponents is clearly larger according to quantum kinetics 
in the first few hundred femtoseconds because of the re- 
tarded onset of the dephasing processes in this theory. 
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10. TR FWM for pulses with \ = 0.3 and delay 
r = 240 fs. B — 18 T. Solid curve: quantum kinetic result; 
Dashed curve: Boltzmann kinetics result. 

However, we point out that the consistence in Fig. 10 
depends on the choice of T and we failed to get consis- 
tent results for all the fields and delays with one constant 
damping T. A more detailed treatment of the damping 
of the memory kernel is needed for a consistent quantum 
kinetic theory. 



Because of the assumed pulse duration of 50 fs is al- 
ready shorter than period of our quantum beats T c — 
2tt/(2lo c ) = 119.7 fs at B — 10 T, one may expect that 
the memory effects of the quantum kinetic scattering in- 
tegrals already may play a role. Therefore we performed 
also quantum kinetic calculation based on Eqs. ( |l2| ) and 
(H). In this case, there is an undetermined damping T 
in the memory kernel. In principle this term should be 
function of fc, q, the magnetic field and the electron dis=. 
tributions, and should be determined self-consistently.ta 
However, the inclusion of such effects will complicate the 
quantum kinetics considerably. Therefore, we took T as 
an adjustable constant. Such a simple approximation 
causes the violation of energy conservation in longtime 



V. CONCLUSION 

In conclusion, we have performed theoretical studies of 
femtosecond kinetics of an optically excited 2D magneto- 
plasmas for intermediate- and high-density excitations. 
Based on a three-subband model in both the conduction 
and the valence band, our study is restricted to high mag- 
netic fields. We calculated the intra- and inter-Landau- 
subband kinetics by bare Coulomb potential scattering 
and found pronounced relaxation oscillations to occur in 
the population of the two lowest Landau subbands for 
lower magnetic fields. We calculated both TR and TI 
FWM signals. Both signals exhibit quantum beats with 
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frequencies around 2lo c . Surprisingly, the resulting de- 
phasing times are rather short and are modulated in the 
strong magnetic field by about 30%. Depending on the 
detuning of the laser pulses, the pulse width, the excited 
densities and the number of excited Landau subbands 
one gets decrease of the dephasing time with increasing 
B field because of the increase of the Coulomb matrix 
elements and the degeneracy of the Landau subbands. 
In regions where the loss of scattering channels exceeds 
these increasing effects, one can also obtain a slight in- 
crease of the dephasing time with the magnetic field. It 
is shown that the retarded onset of the dephasing and re- 
laxation processes in quantum kinetics result in a FWM 
signal which is more coherent in the first few hundred 
femtoseconds. However a consistent theory of the decay 
of the memory kernel in a magneto-plasma is still miss- 
ing. Finally, we point out that more Landau subbands 
are needed to account for the kinetics in lower magnetic 
fields. A corresponding extension of the theory is still 
under investigation and the results will be published else- 
where. 
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Defining two dimensionless variables x = q 2 /4mQ x and 
y = 2xui 2 /n 2 , as Ref. 11, one can express the matrix 
elements in terms of zeroth- (Ko) an d first-order (Ki) 
modified Bessel functions. As there exist the following 
symmetry relations for the matrix elements 



V ■ ■ —V- ■ — (-\\ n+m+i+j V- 



(A3) 



one only needs to calculate 27 terms. The remaining 
terms may be got from the symmetry relations. In the 
following we give these required 27 terms: 

Vfojoo = -e x - y - s K (x) , 
eo 

with 5 = a 2 X 2 (k- k') 2 . 

^io ; oo = --Vye x - y - 3 K (x) . 
eo 

It is noted here and hereafter the term ,/y strictly should 
be y/ysga(q) with sgn standing for the sign function. 
However, due to the fact the all the terms with the factor 
y/y do not appear in the Hartree-Fock terms, we therefore 
just simplify it as y/y. 



Vio ; io = -e x - y -*[(y - x)K (x) + xKi(x)] 
eo 



APPENDIX A: 



Vn;oo = -e x - y -=-[-(x + y)K (x) + xK^x)] 
eo 



We give the analytic expressions for all the Coulomb 
interaction matrix elements Eq. (^|) up to the third Lan- 
dau subbands. We first show the approximation we take 
in our calculation through Voo ; oo- After integrating out 
x and x', Vqo ; oo of Eq. (H) can be written into 



Vr 



00:00 



= ^_ e -±a 2 Sx 2 q - a 2 \ 2 (k-k') 2 

eo 



dz 



Due to the factor e 
increasing \k — k'\ 



— iaX(k 

■\ 2 (k-k') 2 



k')) 2 +q 2 /(2a 2 ) 



(Al) 



Vbo;00 decays quickly with 
Therefore the main contribution to 
the matrix element comes from small values of \k — k'\. 
We therefore neglect the term iaX(k — k') inside the inte- 
grand of Eq. ( Al ) . The integral remaining can be carried 
out analytically. Similar approximation may be applied 
to other matrix elements and therefore, Eq. (^) may be 
written as 



V, 



2we 2 



E 



X <j>n(xi + Sx q )(j)*( x 2 



d{x 1 x 2 )e- iq * {xi - X2) - a2x2{k - k ' )2 

5x q )(f)j(x2)4>m(xi) ■ (A2) 



V w -oi = -e x -y-~[{\ +x- y)Ko(x) - xK^x)] , 
eo 



Vh ; oi = —e x - y - 3 y/y[(l + x + y)K (x) - xK^x)} , 
eo 



Vh iU = — e x - y - 3 {[l + 2x + 2x 2 + y{y -2x- 2)}K (x) 
eo 



+ x(2y - 1 - 2x)K 1 {x)} , 



V 20 ; o = -^ eX v S i( x + y)Ko(x) - xKi(x) 



Vi2;oo — e 

eo 



^[{Zx + y)Ko(x) - 3xKi(x)] , 



^02;0i = -e x -y- s ^/y/2[(y ~ x)K Q (x) + xK^x)] 
eo 
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V 12 -oi = -^e x -y- 3 [(x + y + 2x 2 - y 2 )K (x) - 2x 2 K 1 (x)} , 
v2eo 

e 2 _ 1 1 
V 2;02 = — e x y ^[(-y 2 - xy + x 2 )K (x) + (xy + -x - x^K^x)} , 
e Z Z 

Voo;22 = -e x_w - a [(^ 2 + 3xy + z 2 )^) + - x 2 - 3xy)K 1 (x)] , 
e Z Z 

e 2 1 1 

^20;2i = — e x ^ v ^^/y[(-y 2 + x 2 + x - y - xy)K (x) + (xy - -x - x^K^x)} , 
e a Z Z 

V w ..22 = --e x - y - 3 Vyi(ly 2 -3x-y + xy- 3x 2 )K {x) + (3x 2 - xy + ^K^x)} , 
e Z Z 

v 2 2-,2o = ^- eX ~ v ~ E [(\y 3 - \ x v 2 - \ x2 v + x3 - y 2 + \ x2 + \ x ) K ^ x ) + (J V + \#v - * 3 - \*v - f 2 )^Wl , 

2 /" 

^2o ; oi = -e x - y - H A /|[(y - x - 2)Jf (x) + ^i(x)] , 
e V ^ 

Ki2;io = -4^e x ~ v ~ S [(y 2 -2y-2x- 2x 2 )K (x) + {x + 2x 2 )K 1 (x)] , 
v2e 

e 2 

Vfajii = ^e x " ! '-=[(y 2 - 2xj/ + 2x 2 + 2a; - 2y)aY (x) + (2xy - 2x 2 - x)Kt(x)} , 

v2e 

Vufli = —e x ~ v - 3 M [(3y -3x-2-y 2 + 2xy - 2x 2 )K (x) - (2xy - 2x 2 - 2x)K 1 {x)\ , 
e V 2 

Vu.,22 = --e x -y- 3 {[\{y-2){ y 2 -2 y -2x-xy-2x 2 )+ i -x 2 + 2x*]K {x) + {\xy 2 + x 2 y , 
Eq z z z z z 

V2U22 = i-e°-y-*yfa[{\y* \y 2 + 3xy + \x 2 y \xy 2 + \y - \x - 1 - \x 2 - 2x*)K (x) 
+ [x a + -y x- -x y- -xy + -x + -x )K 1 (x)\ , 

V 02 , 20 = -e x -y- E [(l + 2x -2y-xy + x 2 + ]- y 2 )K {x) + [xy - x 2 - ^K^x)] , 
e Z Z 

V 1 2-20 = --e x -y- E ^{y)[(l + 2x-2y-xy + x 2 + \y 2 )K {) (x) + (xy - x 2 - ^)^(x)] , 
e Z Z 

V 12 ;2i = -e x -y- 3 [(l + 3x-3y- bxy + \y 2 + \x 2 - \y z + 2x 3 - 3x 2 y + lxy 2 )K (x) 
£0 z z z z 

+ + \xy - \x 2 - 2x 3 - ^xy 2 + 3x 2 y)K 1 (x)} , 
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e 2 1 3 7 

V 2 x-2i = — e x - y - a [(-y 3 - 2y 2 + 2y- -xy 2 + Axy - 2x + 3x 2 y - -x 2 - 2x 3 )K {x) 
€q Z Z Z 

+ (^xy 2 - ^xy - 3x 2 y + x + 2x 3 + ^x^K^x)} , 

e 2 1 35 

V 22 -,22 = — e x - y - 3 [(-y 4 - Ax 3 y + 3x 2 y 2 + 2x 4 - xy 3 - 2y 3 + 6xy 2 - llx 2 y + 7x 3 + by 2 - IQxy + — x 2 - Ay 
e 4 4 

9 3 
+ Ax + l)Ko(x) + (xy 3 - 3x 2 y 2 + Ax 3 y - 2x A - -xy 2 + %x 2 y - &x 3 + 6xy - 6x 2 - -x)i^i(x)] . 
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